The electron-phonon theory of superconductors: vertex correction 
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The strong coupling Eliashberg theory plus vertex correction is used to calculate the maps of 
transition temperature (T c ) in parameter-space characterizing superconductivity. Based on these T c 
maps, complex crossovers are found when electron-phonon interaction increases from weak-coupling 
region to strong coupling region. Especially, the interplay interaction between vertex correction and 
Coulomb interaction reduces the higher T c predicted in standard Eliashberg theory to the lower T c 
observed in experiments. The doping-dependent T c of cuprate superconductors can be explained as 
the effect of vertex correction in standard Eliashberg theory and, most importantly, the pseudo-gap 
can be explained as the effect of vertex correction. Our results are also indicated that the non- 
adiabatic effect is the barrier for the realizations of high T c in compounds with compositions of 
light atoms and with high phonon frequencies. 

PACS numbers: 74.20.Fg, 74.25.Dw, 74.62.-c, 71.38.-k 
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I. INTRODUCTION 

In the situation of very strong electron-phonon cou- 
pling, non-adiabatic effects of electron-ion system will be 
so important that the electrons are dressed heavily by 
lattice vibrations and the conventional strong-coupling 
theory [H-[| need be generalized to include the non- 
adiabatic effects or the vertex corrects beyond Migdal's 
theorem. The strong-coupling theory including vertex 
correction has been widely studied by using perturbation 
theory 0-E3| and shown the existences of crossovers when 
electron-phonon interaction evolving from weak-coupling 
region to strong coupling region [IJ, [l2| ■ The problem of 
electron-phonon interaction had been solved numerically 
using Quantum Monte Carlo method (QMC) [l3| and dy- 
namic mean-field theory (DMFT) [li], EH ■ The crossover 
behavior was exhibited in the calculation of Holstein- 
Hubbard model using QMC method [l3| . The crossover 
from electron behavior to polaron behavior with increas- 
ing electron-phonon coupling had also emerged from 
the path- integral calculations in theory of polaron [l6| . 
which had been used to explain the superconductivity 
of copper-oxides high-temperature superconductors by 
forming pairs of polarons p^| - |l8j . 

The standard strong-coupling theory has no bound 
on T c . Recently, the Eliashberg functions a 2 F(uj) ex- 
tracted from the measurements of optical spectroscopy 
for the iron based superconductor Bao.55K .45Fe2As2 [19| 
and the measurements of infrared optical conductivity for 
copper-oxides superconductors [20| predicted very strong 
electron-phonon interaction and very high T c over the 
experimental values [2(j- Additionally, the new found 
rich- hydrogen superconductor silane (S1H4) [2l[ is the re- 
alization of the predicted high-temperature superconduc- 
tor of metal hydrogen by pre-compressed mechanism [22| . 
However, the T c of silane is far lower than the predictions 
of standard string-coupling theory [23| . What is the un- 
derlying reason of so high T c in the predictions of the 
standard strong-coupling theory?. 



In this paper, the T c maps including the influences 
of vertex corrections are studied. Complex crossover 
are found on these T c maps when the parameter A of 
electron-phonon coupling increases from weak-coupling 
region to strong-coupling region. These crossovers are 
very close to the well known A=2 at which the value of 
T c reaches its maximum ||. We also find that T c does 
not monotonously increase with phonon frequency fip. 
T c decreases with fip when Qp higher than a thresh- 
old value. This means that high phonon-frequency is 
unfavorable to superconductivity if vertex correction is 
strong. In section [V] we will show that the vertex cor- 
rection will significantly suppress the predicted T c for 
Bao,55Ko.45Fe2As2 [l9j . the copper-oxides superconduc- 
tors [2(| and superconducting silane 21]. In section ITVl 
main numerical results will be presented. The basic the- 
ory used in this paper will be introduced in section HT1 



II. 



THEORETICAL FORMULAS 



The equation of energy gap in reference [6( had been 
generalized by including the Coulomb interaction [10j. 
The electron-phonon interaction in Nambu's scheme is 
expressed as 



H ep = J2y+T 3 J k ^ k ^ q - k {at qlJ + a q „). (1) 

kqv 



where the index v of phonon mode will be omitted if 
we only consider a single phonon mode. We employ 
the isotropic approximation and the coupling constant 
of electron-phonon interaction is written as J k ,k'v = J- 
The calculations of vertex corrections are greatly simpli- 
fied under isotropic approximation because the electron- 
phonon interactions are included in the vertex corrections 
only by the functions of electron-phonon interaction A„ 
defined as A„ = 2 J °° dcua 2 F(uj)uj/(uj 2 +uj 2 ). The energy- 
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gap equation is expressed as 

\Un\P , 



with the renormalized function Z n ~ Zk(iu> n ), the 
energy-gap function A„ ~ <fik(iuJn) / Zk(iuj n ) and the pa- 



rameters A r , 
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and a„ = (2/-7r) arctan(i?p/.Z„|aj„|). The three parame- 
ters of vertex correction V^ n , , V® n , and C nn > can be found 
in Ref.[l0j]. The Coulomb pseudo-potential is defined as 
M* = Mo/(l + /-iohi(i?s/a;o)), where ^ = N(0)U, U the 
Coulomb interaction between electrons and ujq character- 
istic energy of typical phonon correlated to superconduc- 
tivity. 

If considering that A„ — > when temperature is 
very close to T c , the terms proportional to are 
ignored. The energy-gap equation is simplified to 
Sn'=-oo K nn>{&-n> /\u n >\) = 0. The kernel matrix is ex- 
pressed as 

E-nn' — [A71 — n' E nn ' fl -\- Cnn'\Q"n' ^wi'-^n'j (3) 

TT _ ^n'n"\w n "\ ,\ A 1 

fin' / , I , m 1 Ayi'— n /f -f^n'n" Sn' °ti" Q"n" !• 

n" — — oo 

If the vertex corrections are ignored, three parameters 
Vnn'i Vnn' an d C-nn 1 are all equal to zero and the kernel 
Eq. (|3]) of energy-gap equation reduces to the general form 
without vertex correction Q after some symmetrizations 
and simplifications. In the calculation of a n , Z n ^1 is 
the value of normal state. It's convenient that the K nn > 
matrix is symmetrized as in Ref.Q. The T c is defined 
as the temperature when the maximum of eigenvalues 
-Qmax Q f j cerne i ma trix K nn > crosses zero and changes its 
sign. About iV=200 Matsubara's energies are used to 
solve above equation. Only 20-30 iterations are enough 
to search T c from -600 K to 600 K with accuracy 0.0001 
K by using the bisection method. 

In calculations of a 2 F(u>), the function a 2 (u>) is ap- 
proximately a constant around the peak of phonon mode 
and the density of state of phonon is expressed as 



p^) = <{ ("-fip) 2 +(^) a T^T+WF 1 



- Q P \ < oj 3 
others, 



(4) 



where Sip is the energy of phonon mode, u>2 the half- 
width of peak of phonon mode and W3 = 2u>2- The pa- 
rameter of electron-phonon interaction is defined as 



A = A 



duja 2 F(uj) /uj. 



(5) 



The well known value A = 2 given by McMillan measures 
the instability of superconductivity induced by lattice in- 
stability and plays very important role in this paper. 

Experimental phonon spectrum and the phonon spec- 
trum from linear-response calculation 24] are multi-peak 



structures with broad energy distributions. To make 
the approximation of single peak model to be more re- 
liable approximation, the effective phonon frequency or 
energy should be used. A good choice is the (ui)i n = 
exp(2/A J duj\n(Lj)a 2 F(uj)/u}) defined in Ref. Q. The 
(uj)in is calculated from Eliashberg function that can be 
obtained from linear-response theory or extracted from 
the experimental measurements of phonon properties. 



III. SMALL PARAMETER FOR PERTURBING 
CALCULATION 

It's important to known the small parameter for per- 
turbing calculations in the theory of electron-phonon in- 
teraction. The matrix element of electron-phonon in- 
teraction is given by Jh,Wv = (h/2Alujk-k') 1 ^ 2 ^{k, k') 
(h = 1). The forms of H v (k,k') are dependent on the 
kinds of electron-phonon: deformation potential or polar 
coupling 25]. Under isotropic approximation, the expres- 
sion of electron-phonon coupling constant is given by 



'a si, 



2 N{0) 



(G) 



if the well-known relation MX(lu 2 ) = N(0)(T1 2 ) F m and 
the Einstein spectrum of phonon are used, where Sip is 
phonon energy, A(0) is density of state at Fermi en- 
ergy, Qfm is the averages of Fermi surface and (} is 
the average weighted by (2/X)a 2 F(uj)/uj. The Eliash- 
berg function with Einstein spectrum is expressed as 
a 2 F(uj) — (X/2)ljS(u — Sip). The small parameter for 
perturbing calculation is 



J/E B = 



X Sip 
2h Eb~' 



(7) 



where dimensionless parameter h — N(0)Eb- The stan- 
dard strong-coupling theory is correct only Slp/Ep <C 1 
( Migdal's theorem ). We can see that A > 2 is equiva- 
lent to J 2 > VLpEs/h. The small parameters J/Ep<0-3 
and flp/Es — 0.30 — 0.50 for the fullerides such as 
RkaCgo are estimated by using Eq.([7]) with the param- 
eters A(0)~8 (states/eV), E B ~0.5 cV, A-0.5-1.0 and 
Slp^65-100 meV. Thus strong-coupling theory plus ver- 
tex correct is suitable for fullerides. 

As a comparison, the effective interaction of electrons 
in BCS theory [26j can be written as 



5E( 



2kJg| Jk : k+q\ 2 ^ f t 



kq 



(E k+q - F k y - 



(8) 

Under isotropic approximation, by using uj q ~ Sip, 
Ek+q ~ Ek ~ Ep, v q ~U and Eq.©, the above equation 
can be simplified as 



- Vbcs ' 



,t J 

'k-\-q,a — k — q, 



-k, — aCk,cr 



(9) 



kq 
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where Vbcs = (A — /i*)/2A(0). The small parameter for 
BCS theory is 

Vbcs/Eb = ^f. (10) 

The perturbation theory is correct only (A — n*)/2h < 1 
or A < 2h + fi* , however A > fi* to keep supercon- 
ducting state stable. If A > 2h + /i*, the strongly 
coupling electron-pairs will significantly modify the elec- 
tronic structure of superconductors and lead to structural 
instability. 

IV. THE GENERAL RESULTS OF VERTEX 
CORRECTIONS 

The three-dimensional Tp maps had been calculated 
in the previous paper [l0|. In this section, three- 
dimensional T c maps in X-flp-fi* phase space includ- 
ing vertex corrections are calculated by using the simple 
phonon spectrum with the form of Eq.Q. From electron 
point of view, the vertex correction or non-adiabatic ef- 
fect can be controlled by the effective width E^, on the 
other hand, from ion point of view, it can be controlled 
by the cutoff of phonon energy ujq. In this work, the ver- 
tex corrections are controlled by the effective band-width 
Eb with the range from 0.5 eV to 5 eV. The situation 
Eg = oo is equivalent to no vertex correction. 

The FigfIJa,b,c) illustrate the evolution of T c map with 
decreasing Eg. The FigfTJa) is the T c map having been 
obtained in the previous work without considering vertex 
corrections [Io| . The large deformation with strong ver- 
tex correction starts at E#=1.7 eV in Figfljb) near the 
well known A=2.0 and high phonon energy. With E# 
decreasing to 1 eV further, the region with strong ver- 
tex correction rapidly expands and occupies large part 
of parameter space with fip>80 meV in FigfTJc). In 
the region with J7p<80 meV, the T c is strongly sup- 
pressed however there are no discontinuous changes of T c 
or breaking of contour lines. An important result from 
the FigfTJis that T c does not change with A monotonously 
if phonon energy ftp is high enough. The FigJ^a) shows 
the changes of T c with A along two arrows (A,B) shown 
in FigfTJc). If ttp=80 meV, the T c monotonously in- 
creases with A. However for fip=90 meV, the T c first 
increases with A, reaches the maximum at A ~ 1.5 — 1.7 
and then quickly decreases with increasing A. Further 
increasing A>2, T c will be very low due to strong ver- 
tex corrections. The non-monotonous A-dependent T c in 
Fig|5][a) had been found in the non-adiabatic theory of 
superconductivity [111 ]. Some crossover behaviors from 
weak coupling to strong coupling region had been pre- 
dicted in Holstcin-Hubbard model solved numerically by 
quantum Monte Carlo method (l3| and in polaron the- 
ory [l6j]. It's very reasonable that the non- monotonous 
A-dependent T c is equivalent to the crossovers found in 
QMC calculation [l3[ and polaron theory [l6j]. So only 
the leading vertex correction can describe qualitatively 
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FIG. 1: The evolution of T c map on A-Slp plane with in- 
creasing strengths of vertex corrections (decreasing effective 
band- width Eb) with (a) Eb = oo, (b) Eb = 1.7 eV and (c) 
Eb=1 eV. The Coulomb pseudo-potential fi* =0.1. 
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FIG. 2: (a) The T c change along two arrows shown in FigtTJc) 
with fixed phonon energies S7p=80 meV and 95 meV respec- 
tively, (b) The T c change along two arrows shown in Fig[3jc) 
with fixed Coulomb pseudo-potentials /x*=0.10 and 0.25 re- 
spectively. 
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FIG. 3: The evolutions of T c map on p*-A plane (fip=69 
meV) with decreasing effective band-width (a) E_g = oo, (b) 
E s = 1.7 eV and (c) E B = 1.0 eV. 



very well the electron-phonon interaction in strong cou- 
pling region. 

In the previous paper, the Coulomb pseudo-potential 
H* had small effects whenu*>0.2 in the calculations 
without vertex corrections [lfj]. The Fig|3][a) is the nor- 
mal T c map on /j,*-X plane without vertex correction [To| . 
The figure shows that when /i*>0.2, T c is insensitive to 
the change of /i*. The breaking contour line with T c =0 K 
are because of the inaccurate calculations when T c <0.1 K 
if only N—200 Matsubara energies are used. The contour 
lines with T c >0.1 K is accurate enough. If the Coulomb 
pseudo-potential and vertex correction work together, 
the situation will change drastically and some new in- 
teresting results will appear. The large deformations are 
found in FigfSJc) if decreases to 1.0 eV. In the re- 
gion with A<2.0, the contour lines with higher T c are 
moved out but the contour lines with lower T c fill in this 
region. As expected, the large deformations and discon- 
tinuous changes of contour lines appear on the T c map 
when /i*>0.20. The contour lines with iso- values from 
T c =20 K to 200 K are packed together within the region 
0.20<^*<0.25 and A>2. The figure clearly shows that 
if the Coulomb pseudo-potential is larger enough, the 
T c will change with A non-monotonously. The changes 
of T c along two arrows with /^*=0.1 and 0.25 are plot- 
ted in Figl^b). For ^*=0.25, T c first increases with A 
until reaches the maximum at A=2.2 and then sharply 
decreases to smaller value at A=2.5. The crossover be- 
havior is enhanced by strong Coulomb interaction. 

There are three methods that can control vertex cor- 
rection: (1) only wo changes and E# keeps unchanged, 
(2) Eb changes and wo keeps unchanged , (3) both Ep 
and wo change. It's possible that three methods provide 
qualitatively different results. In this paper, we use nar- 
row peak approximation so that wo ~ flp. The T c map 
on fl P -E B plane is plotted in Figg^b). The FiggJJa) 
shows how T c changes with Vlp/Ep along the direction 
of arrow A in FigQJb) with Es=l eV. It's found that at 
first T c increases with Vtp/Ep and then decreases with 
it. The arrow A in Fig01Jb) is corresponding to the hor- 
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FIG. 4: (b) T c map on E B -Slp plane with //=0.1 and A=1.0. 
(a) T c changes non-monotonously with increasing Qp/Ep 
along direction of arrow A with increasing tip and fixed 
E B = 1.0 eV. 



izontal arrow C in FigQJc). The behavior of T.-with 
ujq/Ef had been found in a series of references @,H|, EI 
in that coupling constant of electron-phonon interaction 
9k,k+q ( or Jk,k+q m this paper) is dependent on a cut- 
off q c for wave- vector q. From FigQJc) an d Fig|4jb), the 
increasing T c at small flp/Ep isn't the effect of vertex 
correction and is the standard result of strong-coupling 
theory. The real effects of vertex corrections are that T c 
decreases with VLp/Ep at large Vtp/Ep. Certainly, the 
vertex correction is enhanced if Qp / Ep increases by de- 
creasing E^ with fixed tip, illustrated by the horizontal 
arrow B in Fig[4jb). 

Finally, T c map on E^-A plane is presented in Fig[S] 
with flp— 72 meV. If E^ increases but A keeps un- 
changed, the T c monotonously increases with E# until 
to the limit of non- vertex correction. More interestingly, 
on this map, the T c is non- monotonous dependent onEg 
along strait line from PI to P2 companying by the de- 
crease of A from 3.0 to 0.2. The non-monotonous depen- 
dence of T c on effective band- width E^ is similar to the 
band-filling effects of T c in the non-adiabatic theory of 
superconductivity [9]. Our results show that, if f2p >80 
meV, the suppression of Tc will be more prominently. 



V. DISCUSSION 

As mentioned in the introduction, the values of T c ob- 
tained from standard strong-coupling theory are gener- 
ally higher than those measured in experiments. The 
copper-oxides superconductors Bi2Sr2CaCu20s+5 and 
Bi2Sr2Cu206+5 studied in Ref 20] have very strong 
electron-phonon interactions A^2.36-2.85 and overesti- 
mated T c in underdoped samples. With increasing dop- 
ing 5, the values of A decrease to about 0.35-1.42. The ef- 
fective band- widths Es of conducting electrons for these 
cuprates are distributed from 1 eV to 3 eV. The effective 
phonon energies are distributed from 50 meV to 80 meV. 
The Coulomb interactions are strong in cuprate super- 
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FIG. 5: (a) The T c map on Eb-X plane with /n*=0.25 and 
Qp=72 meV. (b) The open circle line is the evolution of T c 
from PI to P2 in (a). The solid line T* is the standard results 
in strong coupling theory without vertex correction, (c) The 
5-\ relation is adopted in Ref. 20 



conductors fi*—0.25. We can clearly see from FiglJfb) 
and Figj3fc) that in underdoped samples, T c can be re- 
duced to very small values in strong clcctron-phonon in- 
teraction region (A>2.5) due to the interplay interaction 
between vertex correction and strong Coulomb interac- 
tion. As shown in FigJSJb), the values of T c for all sam- 
ples are reduced from around 200 K to lower than 150 
K and close to experimental values [2(J. Our results are 
provided an explanation to pseudo-gap in underdoped re- 
gion shown in Fig|5]Jb). The cooper-pairs pre- form at the 
mean field value T* of transition temperature in standard 
strong-coupling theory. However strong non-adiabatic ef- 
fects induce the instability of cooper-pairs and suppress 
real T c to lower values. The T* degenerates with T c in 
overdoped region is similar to the example (1) of Fig. 9 in 

For the iron-based superconductors Bao.55Ko.45Fe2As2 
studied in Ref. fl9| , the parameters of boson-phonon in- 
teraction decrease from 3.42 to 0.78 with increasing tem- 
perature. The contributions of electron-phonon interac- 
tions have energies lower than 40 meV in Eliashberg func- 
tions a 2 F(uj). At 28 K close to T c , the roughly estimated 
value of A contributed by energies lower than 40 meV is 
about 2.0 close to the values at temperature 86 K. This 



is too strong to account to experimental T c =28 K. Be- 
cause the effective phonon energy is about 20-40 meV, 
the reduced T c from the vertex corrections are also small 
(Fig|T} although the effective band-widths of conducting 
electrons are small from 1 cV to 2 eV. The key problem 
is that how we can extract the more accurate value of 
A of electron-phonon interaction from total A. If we use 
the A=1.44 obtained at temperature 151 K and the effec- 
tive phonon energy Jlp=20 meV, the obtained T c =29 K 
is close to experimental T c =28 K. The reasonable value 
for iron-based superconductor is about A=1.0. This is 
indicative that the contribution of electron-phonon in- 
teraction to Eliashberg function a 2 F(uj) dominates over 
other interactions at high temperature. 

From Fig|4] we can see that if effective band-width 
Es is small, the T c will decrease with increasing phonon 
frequency if Qp is higher enough. This means that 
we can not increase T c infinitely by decreasing mass of 
atom or increasing phonon frequency. The strong ver- 
tex corrections in high phonon-energy region suppress 
T c to smaller values than those of predicted by standard 
strong-coupling theory. This is a reasonable explanation 
for the lower T c in the experiment of silane superconduc- 
tor [2l[ than that in the strong-coupling theory (23|. 

The superconductivity of fullerides can be described 
using Cooper-pairs glued by virtual phonon excitation. 
However there are some very important problems are 
unsolved (28j . The fullerides have very large number 
of phonon modes from low frequency to high frequency. 
What phonon modes are most important to their super- 
conductivity?. There are two groups of phonon modes 
of fullerides: (a) inter-molecular molecular modes with 
energies lower than 150cm -1 (or 18 meV) and (b) intra- 
molecular phonon modes with energies higher than 250 
cm -1 (or 31 meV). The inter- molecular modes have small 
contributions to T c due to their lower energies and not 
too strong electron-phonon interaction. Most impor- 
tantly, our results shown in FigOJc) and Fig0]prove that 
the 'tangential' intra-molecular modes with energy higher 
than 1000 cm -1 (or 124 meV) have small contributions 
to T c because of strong vertex corrections. So the intra- 
molecular 'radial' modes from 250cm -1 to 1000cm -1 of 
group (b) should contribute main parts of superconduc- 
tivity. In fact effective phonon frequency (u>)i n is located 
in this energy range. The empty ti M orbits will form nar- 
row energy bands with width about 500 meV in Cqq solid 
and strongly couple with intra-molecular H g phonons. In 
RbaCgo solid, the width of ti u bands increase to about 1 
eV and is half filled with electrons. Based on our calcu- 
lations using density functional theory with plane-wave 
pseudo-potential methods, the effective band- width Eb is 
higher than 500 meV and lower than 1 eV. If we choose 
the intra-molecular radial mode f2p=525 cm -1 (or 65 
meV) and the half- width Ep=500 meV, T c will be lower 
than experimental value due to vertex-correction effects 
shown in Fig[^a). If we slightly increase Eb and use the 
averaged value Ep=750 meV, both T c =29.5 K and iso- 
tope coefficient a ~ 0.30 [29[ are close to the values in 
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experiments if A=1.04 and fx*=0.1. The required param- 
eter A is different from our previous work value 0.66 [3(| 
in the calculation without vertex correction, however still 
close to the reasonable range from 0.5 to 1.0 [28]. The 
accurate calculations of A and /i* are important to un- 
derstand the role of vertex correction. 



VI. CONCLUSION 

In summary, the strong-coupling Eliashberg theory 
including vertex correction is systemically studied in 
this paper. The T c maps in parameter-space X-fl-fi* 
for different contain the completed information on 
strong-coupling theory and the vertex corrections. Es- 
pecially, the combined interaction of vertex correction 
and Coulomb interaction can significantly depress T c 
to small value. The non-monotonous changes of T c 
with increasing A show the crossover behaviors near A=2 
when A evolving from weak-coupling region to strong- 



coupling region. The crossovers can explain the doping- 
dependent T c of cuprate superconductors if the Coulomb 
interactions are strong. The crossover behavior in In- 
dependent T c indicates that high-frequency phonon is 
unfavorable to high T c . Thus the strong non-adiabatic 
effect makes it hard to find high T c superconductors in 
compounds containing light elements. Finally the T c 
maps in the previous paper [10( and the maps with ver- 
tex corrections in this paper provide very comprehensive 
understanding of superconductivity of superconductors. 
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